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Stretched exponential relaxation of a quantity n versus time t according to n = no exp[— (A*f)''] 
is ubiquitous in many research fields, where A* is a characteristic relaxation rate and the stretching 
exponent /3 is in the range < /3 < f . Here we consider systems in which the stretched exponential 
relaxation arises from the global relaxation of a system containing independently exponentially 
relaxing species with a probability distribution P(A/A*,/3) of relaxation rates A. We study the 
properties of P(A/A*,/3) and their dependence on 13. Physical interpretations of A* and /3, derived 
from consideration of P(A/A*,/3) and its moments, are discussed. 
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INTRODUCTION 

The stretched exponential relaxation function de- 
scribes the time t dependence of a relaxing quantity n, 
according to 



n = no exp[— (A*i)''] , 



(1) 



where no = n(t = 0), A* is a characteristic relax- 
ation rate, and the stretching exponent (3 is in the range 
< /3 < 1. This behavior has been observed for a wide 
variety of physical quantities in many different systems 
and research areas. |l| Log-linear and linear-log plots of 
the stretched exponential function versus \*t, which are 
the two most common ways of plotting the stretched ex- 
ponential function, are shown in Figs.^a) and^b), re- 
spectively, for (3 values from 0.1 to 1 in 0.1 increments. 
All of the plots cross at a time ti/g = 1/A* at which the 
stretched exponential function has the value for all (3. 
A pure exponential decay, corresponding to /? = 1, plots 
as a straight line in Fig.^a). As (3 decreases from unity, 
more and more pronounced positive curvature is evident 
at small t. At small times, a Taylor series expansion of 
Eq. Q for A*t <C 1 gives 
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no 
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(2) 



Thus the stretched exponential function with < /? < 1 
is singular at t ~ 0, with an infinitely negative slope 
there. 

A natural and commonly used interpretation of an ob- 
served stretched exponential relaxation is in terms of the 
global relaxation of a system containing many indepen- 
dently relaxing species, each of which decays exponen- 
tially in time with a specific fixed relaxation rate A. Then 
one can write the stretched exponential function as a con- 
tinuous sum of pure exponential decays, with a particular 
probability distribution P of A values for a given value of 
f3. In such a probability distribution, one must normal- 
ize A to the characteristic relaxation rate A* appearing 
in the stretched exponential function Hence for such 
systems one can write the stretched exponential function 
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FIG. 1: Log- linear (a) and linear- log (b) plots of the stretched 
exponential function in Eq. Q versus X*t for /3 values from 
0.1 to 1 in 0.1 increments. 



in Eq. Q as 



where 



(3) 



Since P{s, P) is a probability density, one has 
/„°°P(s,/3)ds = l. 
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A stretched exponential form of the ''Li nuclear 
spin-lattice relaxation following saturation was recently 
observed[3, iii ^Li NMR experiments on the 

heavy fermioiij5( compound LiV2 04 containing mag- 
netic defects, y] In our attempts to understand the ori- 
gin of this nonexponential relaxation, we considered the 
above model of independently relaxing nuclear spins and 
found the published information on P{s,(3) (e.g., Refs. 
HSBHEIilElEltobe insufficient for our pur- 
poses, particularly with regard to the systematic behav- 
iors of P{s) versus P and to the physical significances 
of the parameters A* and /3 in the stretched exponen- 
tial function Since in the case of independently re- 
laxing species the essential physics of the system resides 
in P{s, f3) rather than in the relaxation function it- 
self, we considered it important to further study P(s,/3). 
Here we report the results, some of which were briefly 
mentioned in Ref. 0. 



RESULTS 
Probability density 

For /9 = 1, the stretched exponential function is a 
pure exponential with relaxation rate A = A* and hence 
the probability density P{s, 1) in Eq. is a Dirac 6 func- 
tion at s = 1. For general /?, from Eq. Q one sees that 
P(s, f3) is the inverse Laplace transform of the stretched 
exponential, given by 



(4) 



The change of variables u = —ix allows one to write 
P(s, (3) as the Fourier transform 
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One can also express P{s,(3) as|6l ISllllI 



sin(n7r/3) 



(5) 



(6) 



where T{z) is the Gamma (factorial) function. Using 
Mathematica 4.0, we have obtained from Eqs. lO and 
© closed analytic solutions for P{s, (3) for many rational 
values of f3. The simplest expressions are obtained for 
f3 — 1/3, 1/2, and 2/3, for which solutions have also 
been given in Refs. andlioi For /? = 1/3, we find the 
two alternative expressions 



P{s,-) = iz^M{z) 



37rs'^/ 



(7) 



where Ai(z) is the Airy function with z = (3s) and 
Kn{x) is the modified Bessel function of the second kind. 
For 13 = 1/2, we get 



/In 1 /In 

Pis, -) = , exp ( - — ) , 



(8) 



which is also listed in tables of Laplace transforms. |l4| 
For [3 = 2/3, we obtain 

P{s, ^) = 6z^/4 exp ( - ^) [Ai(z) - ^ Ai'(z)] , (9) 

where z = (3s) 

In general, the analytic results for rational j3 = n/m 
values are expressed by Mathematica in terms of general- 
ized hypergeometric functions pPg(a; b; z). For example, 
for (3 = 2/3, P(s, (3) can also be expressed as 



P 



(4) 
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(10) 



with z = — 4/(27s^). The expressions become progres- 
sively longer as the integer denominator m increases. 
For example, the expression for P{s,(3 — 4/5) contains 
four additive pPg(a;b;z) terms with p = q ~ 3 and 
z = -256/(31255-^), with l/s'^+"''^ multiplicative pref- 
actors where n — 1-4, respectively. 

The probability density P{s,P) in Eq. ^ is real by 
definition. By finding the real and imaginary parts of 
the integral on the right side of Eq. , one can explicitly 
show that the imaginary part is identically zero. The real 
part, which is P(s, /?), is found to be 

1 f°° 

TT Jo 

This integral can be evaluated numerically for arbitrary 
(3 using Mathematica. Similar expressions were obtained 
in Refs. i and O 

The evolution of P(s, (3) versus s at fixed (3 for several 
rational values of (3, from our analytic results, is plotted 
in Fig. Ela) on linear scales and in Fig. Etb) on semilog 
scales. The same type of plots as in Fig. EJa) were pre- 
viously given in Fig. 1 of Ref. IT3. and related plots were 
given m Refs. H andE From Fig. Efa), a (3 < 1 causes 
the infinitely high and narrow Dirac 5 function probabil- 
ity distribution for /3 = 1 at s = 1 to broaden. P(s) be- 
comes highly asymmetric, and the peak in P(s) becomes 
finite and moves towards slower rates which is compen- 
sated by a long tail to faster rates. 

The value of s at which P(s,(3) is maximum for fixed 
/3, s{P'^^^), is plotted versus /3 on linear and semilog 
scales in Figs. IH^a) and (b), respectively. We find that 
s(P™^'^) decreases with decreasing f3 and approaches zero 
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FIG. 2: Linear (a) and semilog (b) plots of the probability 
density P(s,/3) for the relaxation rate in Eq. ^ versus nor- 
malized relaxation rate s = A/ A* for several rational values 
oip. 




FIG. 3: Linear (a) and semilog (b) plots of the value s(P™="') 
of s at which the probability density P(s) in Eq. Q (see Fig. 
1^ is maximum, plotted versus the stretching exponent /3. The 
solid curve in (a) and (b) is an exponential fit to the data for 
P < 0.5 as given in Eq. 1121 . The position of the peak can be 
taken to be the intrinsic small-s cutoff to P{s) (see text). 



exponentially for f3 < 0.5. A fit to the s{P™'^^) versus f3 
data in the range l/8</3<l/2 yielded 

^(^'-)-P-P(-^)- (12) 

Plots of this expression are shown as the solid curves in 
Figs. Eta) andEtb). 

For small s, there is an exponential decrease in P{s, (3) 
with decreasing s that is given in the limit s — > (i.e., 
below the peak) by|l ITU 

P{s ^ 0, /?) = exp ( - ^^7(^) , (13) 

where 

/3l+/3/[2(l-/3)] 
yW(l^ 

and 

6 = (1 - /3)/3Mi-« . 



From Eq. for /3 = 1/2 one gets directly that a ~ 
l/-\/47r and h = 1/4, in agreement with these results. 
Thus the probability distribution P{s,(3) contains an 
intrinsic low-s cutoff that decreases with decreasing f3. 
Note that the peak position of P{s, (3) in Eq. (|12() has 
the same form as the s — > behavior in Eq. 113|) . We 
take the value of the low-s cutoff to be the position of the 
peak in P(s, (3) that is plotted in Fig.|3|and described by 
Eq. ini) for /? < 1/2. 

In the opposite limit of large s 1, in Eq. lO one 
retains only the first term in the power series, giving 

P(s^oo,/3) = ^ , (14) 

where 

c = £i/?±ll sin(7r/3) . (15) 

TT 

For < /? < 1, one has that r(/3 + 1) - 1. 

The median Smedian of the probability distribution 
P{s) has not been discussed before in the literature to 
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FIG. 4: (a) Median of the probability distribution P{s) versus 
13. (b) Expanded plot of the data in (a) for 0.3 < /? < 1. The 
lines are guides to the eye. 
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FIG. 5: The integral P(s,/3)ds versus 



our knowledge. It is defined to be the value of s for 
which it is equally likely for s to be less than Smodian 
as it is to be greater. It is calculated by solving the 
expression J^"""^"^'' P{s, f])ds — 1/2 for a fixed value of 
f3. The Smcdian is plotted versus /3 in Fig. ^a), and an 




FIG. 6: Probability that s is less than the value so, given by 
the integral J^" P(s,/3)ds, for /? = 1/8 and 7/8. 



expanded plot of the data for 0.3 < /? < 1 is shown 
in Fig. EJb). One sees from Fig. 0] that with decreas- 
ing /3, Smedian remains nearly equal to unity from 13 — 1 
down to about j3 = 0.5, below which Smodian begins to 
increase dramatically. A related quantity is the integral 
P{s, j3)ds^ which measures the probability that s is 
less than or equal to unity. If the median is at s = 1, 
then the integral should equal 1/2. The integral is plot- 
ted versus (3 in Fig. El In the range 0.1 < /3 < 0.9, the 
integral is equal to 1/2 to within ±0.1. The reason for the 
difference between the (i dependences of the two quanti- 
ties is apparent from plots for (3—1/^ and 7/8 in Fig. |H| 
of the probability that s is less than a value sq, given by 
the integral J^" P{s,f3)ds. One sees from the figure how 
it happens that for /3 = 1/8, Smcdian > 2 P{s, P)ds, 
whereas for (3 = 7/8, Smcdian < 2 P{s,l3)ds. 



Moments of P{s,l3) 

The moments of the probability distribution are de- 
fined by (s")avc = s"P{s,(3)ds. From Eq. (dj), the 
n = 1 moment, which is the average Savc = (A/A*)avc — 
Aavc/A*, is infinite for all 13 with < /? < 1. Hence the 
average relaxation rate Aavo is infinite. Indeed all positive 
moments such as also (s^)avo are infinite. Very fast (infi- 
nite) relaxation rate components are required to produce 
the infinitely negative slope for any /3 < 1 in the stretched 
exponential function (QJ as t ^ as described in Eq. (0) 
and seen in Fig. Q^a). A cutoff to P(s,/3) at large s is 
required to obtain a finite initial slope of n{t) in Eq. 
and finite averages (s")ave- Every real system must have 
a high-s cutoff to P(s, but of course it will depend on 
the system under consideration. We will return to this 
issue in Sec. g 

On the other hand, the n = —1 moment of P{s,f3), 
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FIG. 7: Semilog plot of the average (A*/A)avc versus the 
stretching exponent 13 according to Eq. H16|l with m = 1. 



which is the average (l/s)avc — (A*/A)avo, is finite for 
< /? < 1 and increases monotonically and rapidly with 
decreasing /3, diverging at /3 = 0. In general, one hasQ 



T{m/(3) 



(16) 



where m = —n. A semilog plot of (A*/A)avc versus [3 for 
0.01 < /3 < 1 is shown in Fig. 



Physical interpretations of A* and (3 

The parameter A* in the stretched exponential function 
is often referred to in the literature as some undefined 
"characteristic relaxation rate" or as some undefined "av- 
erage" relaxation rate. As shown above, A* is neither the 
average of A [which is infinite in the absence of a high-s 
cutoff to P(s,/3)] nor the inverse of the average of 1/A. 
Evidently A* is a characteristic property of P{s,/3) itself 
and not of its moments. Indeed, from the above discus- 
sion and the data in Fig. we infer that the physical 
interpretation of A* is that A is about equally likely to be 
less than A* as it is to be greater (to within ~ ±20%). 
From Fig. 21 one sees that A* is within about 10% of the 
median of P(s,/3) for 0.5 < /3 < 1. For /3 < 0.5, the 
median strongly increases due to the long high-s tail to 
P{s,f3), but the integral in Fig. continues to decrease 
slowly with decreasing /3. 

The stretching exponent (3 is often cited as a measure 
of the width As of the distribution P{s,f3). However, a 
statistical definition of the width such as the rms width 
As = (s2)avc — (save)^ IS Undefined for P{s,(3) since 
as shown above both of the averages in the square root 
are infinite in the absence of a high-s cutoff to P(s,/3). 
Another possibility is that /3 is a measure of the full width 
at half maximum (As = FWHM) of P(s,/3). From Fig. 
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FIG. 8: Full width at half maximum peak value (FWHM) 
versus f3 of the probability distribution P(s) on a linear s 
scale (a) and of P(logjQ(s)) on a logarithmic s scale (b). See 
Figs. |2Ia) and (b), respectively. The solid curves are guides 
to the eye. 



IHl^a), the FWHM initially increases as (3 decreases below 
unity, but then decreases as /3 decreases further for /3 < 
0.65. Therefore, /3 is a multivalued function of FWHM, 
and hence correlation with the FWHM is not a useful 
physical interpretation of f3. On a logarithmic s scale, 
the FWHM increases monotonically with decreasing /3 
as shown in Fig. |Htb), and appears to diverge as /3 — > 
0. Therefore /? can be considered to be related to the 
logarithmic FWHM of P(s), but not to the FWHM itself. 

Therefore we seek a physical interpretation of (3 that 
is not exphcitly tied to the width of P(s, /3) versus s. We 
have seen above that the position s{P'^^^) of the peak 
in P{s,(3), plotted in Fig. |31 versus /3, decreases mono- 
tonically with decreasing and that below the peak, 
P(s, /3) decreases exponentially as s decreases. As dis- 
cussed above, this demonstrates that P{s,(3) has a low-s 
cutoff that decreases monotonically with decreasing (3. 
Thus a useful physical interpretation of (3 is that (3 is 
a measure of the intrinsic small-relaxation-rate cutoff of 
P{s,/3). To our knowledge, this identification and the 
above physical interpretation of A* here and in Ref. 
have not appeared elsewhere in the literature. 
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One can rewrite the stretched exponential function 



no exp[-(i/r* 



(17) 



where the relaxation time of a particular relaxing species 
in the system is r = 1/ A and the characteristic relax- 
ation time of the probability distribution is r* = 1/A*, 
so r/r* = A*/A. From the above results, one obtains, 
for example, that the plot of (A*/A)avc versus (3 in Fig. 
[7|is the same as a plot of Tavo/T* versus /?. From Fig. 
the normalized average relaxation time t.^vc/'t* is well- 
defined at each /3, increases monotonically and rapidly as 
(i decreases below unity, and diverges as ^ 0. Thus 
an additional physical interpretation of (3 is that 1/(3 is 
a measure of the average relaxation time of the relaxing 
species in the system relative to the parameter r* = 1/ A* . 



Influences of a large-s cutoff to P{s,l3) 

We have alluded above to the singular nature of the 
stretched exponential function at time t = 0. Asso- 
ciated with this singularity are infinite averages (s")avG 
with positive integer n for any fixed P with < /3 < 1. 
Here we briefly discuss the influences of a large-s cutoff 
■Scutoff to P(s, /?) on Save ^ud on the relaxation function 
that can be computed from the modifled P{s, 13, Scutos)- 
That relaxation function is no longer a pure stretched 
exponential. By deflnition, then. 



Pis, 13, Scutoff) =0 for s > s, 



cutoff 



(18) 



For illustrative purposes, we consider the probability 
distribution for (3 — \/2 because of the simple form of 
P{s,(3 = 1/2) in Eq. ((SJ. This (3 value is also commonly 
observed in real systems so the discussion here may have 
practical applications. The only influence of the cutoff on 
P(s, f3), beyond the cutoff condition ((TH|) . is that P(s, (3) 
must be renormalized so that 



P(s,/3, Scutoff )rfs = 1 



(19) 



Applying conditions and to P{s, 1/2) in Eq. 
gives, for s < Scutoff, 



Pis, 1/2, Scu 



toff) 



exp(-i7) 



vItts^ erfc( ., , ^ ) 



(20) 



where erfc(2;) is the complementary error function. 

The average of s is Savc = Jq""'''''* sPis, f3, Scutos)ds, 
yielding for /3 = 1/2 



V^cutoff exp ( 



4s<,„toff ' 



/tt erfc( 



Zy'Scutott 



) 



(21) 
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FIG. 9: Linear (a) and log-log (b) plots of Save = Aavc/A* 
in Eq. 12H versus the large-s cutoff Scutoff to the probability 
distribution P(s,/3 = l/2,Scutoff) in Eq. iCTl . 



Linear and log-log plots of Savo versus Scutoff are shown 
in Figs.inia) andEl^b), respectively. For Scutoff ^ 1, one 
can approximate Eq. H21|) as 



1 1 



-^cutoff 
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(Scutoff > 1) 



From Eq. Q , one obtains the time dependence of the 
relaxation of the quantity n(t) using 



nit) 

no JO 
or, with a cutoff present, 

nit) 



P{s,P)e-''^'*ds 



no 



/ P(s,/3, Scutoff )e" 



sA*t 



ds 



(22) 



(23) 



Shown in Fig. llUf a) are plots oin/no versus t for (3 — \/2 
and for Scutoff = 100 and oo, where P(s, 1/2, Scutoff) is 
given in Eq. (|20|l . Expanded plots for A*t <^ 1 are shown 
in Fig. lior b). If there is no large-s cutoff, then Savc = oo 
and the initial slope oin/no versus t is — oo for < /? < 1 
as previously discussed. However, the presence of the 
cutoff gives a finite Savc and an initial linear decrease of 
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FIG. 10: (a) Comparison of the time t dependent stretched 
exponential relaxation with (3—1/2 without a cutoff to the 
probability distribution and with a cutoff Scutoft = 100. (b) 
Expanded plot of the data in (a) near t = 0. Note the initial 
linear dependence of n on t with the finite cutoff Scutoff = 100. 



n{t)/no versus time as seen in Fig. [TUT b'). In that case, 
at sufficiently small t « (scutoff A*)^^ one can Taylor 



-sX't 



1 - 



expand the exponential in Eq. (|23|l to get e 
sX*t. Substituting this expression for the exponential 
into Eq. gives 



n{t 0) 



3A*t — 1 — Aavc^ 



no 



where we have used Save — Aavc/A*. The initial slope 
of n/no versus t is then a true measure of the negative 
of the average relaxation rate Aave m the system. Every 
real system must have a large-s cutoff. Experimentally, 
an important resolution issue is being able to look at 
short enough times to be confident that one is measuring 
the initial slope of the relaxation. The length of time 
that the initial linear relaxation is retained increases as 
Scutoff decreases. 



SUMMARY 

The ubiquitous stretched exponential relaxation func- 
tion has been found to apply to the relaxation behav- 
ior of many different systems. In this paper we have ex- 
amined some systematics of the probability distribution 
P(s,/?), where s = A/A*, that apply in the particular 
case that the stretched exponential function arises from 
the global sum of exponential decays of independently 
relaxing species with relaxation rates A. 

The functional dependence of the peak position of P{s) 
on P has been determined. This peak position decreases 
monotonically with decreasing /3 and was characterized 
as a measure of the intrinsic low-s cutoff possessed by 
P(s). We therefore suggested that a physical interpre- 
tation of P is that /3 is a measure of this intrinsic low-s 
cutoff. Additionally, 1//3 is a measure of the average re- 
laxation time (not rate) of the relaxing constituents of 
the system. These interpretations contrast with a com- 
mon one that /3 is a measure of the width in s of P(s, 
which was shown to be not quantitatively useful. 

We derived and discussed the /3-dependence of the me- 
dian of P(s). The positive moments of P(s), which are 
the averages (s")avo, are infinite in the absence of a large- 
s cutoff to P(s, /?), so A* is not related to any such average 
as often assumed. We suggested instead that the funda- 
mental physical interpretation of A* is that A is about 
equally likely (to within ~ ±20%) to be less than A* as 
it is to be greater. 

The influence of a large-s cutoff to P{s,(3) on Save = 
Aavc/A* and on the time dependence of the relaxation 
were investigated for the illustrative case of /3 = 1/2. 
The cutoff leads to a finite average relaxation rate that 
decreases as the cutoff decreases. The cutoff also removes 
the infinite-slope singularity in the stretched exponential 
relaxation at time t = and replaces it with an initial 
linear time dependence, the slope of which is the negative 
of the average relaxation rate of the constituent relaxing 
species in the system. 

We note that experimental time-dependent relaxation 
data can often be fitted by a discrete sum of exponen- 
tial relaxations, as well as by the stretched exponential 
function. Thus experimental resolution is important to 
distinguishing between such similar types of fits. Ulti- 
mately, the justification for using the stretched exponen- 
tial function to fit relaxation data, as opposed to some 
other function such as a discrete sum of decaying expo- 
nentials, must reside in the physics of the system under 
consideration. 

The author is grateful for the collaboration of and 
discussions with his coauthors of Ref. 4. Ames Labo- 
ratory is operated for the U.S. Department of Energy by 
Iowa State University under Contract No. W-74G5-Eng- 
82. This work was supported by the Director for Energy 
Research, Office of Basic Energy Sciences. 
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